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EXECUTIVE  SUMMARY 


This  report  defines  a  general  procedure  for  the  reduction  of  center  of  gravity 
vertical  acceleration  data  towards  gust  statistics.  The  specific  aspects  treated 
include: (1)  Separation  of  accelerations  due  to  gusts  and  maneuvers. 

(2)  Definition  of  peaks  and  valleys. (3)  Reduction  of  accelerations  to  "gust 
velocities",  using  either  a  discrete  or  a  continuous  gust  concept. 

The  possibility  of  reevaluating  and  reducing  existing  European  Load  Data  sources 
into  a  format  compatible  with  the  defined  general  reduction  procedure  is 
investigated . 


L. _ INTRODUCTION 


As  part  of  the  research  program  on  aging  aircraft,  the  Federal  Aviation 
Administration  (FAA)  has  started  an  ambitious  program  on  flight  loads.  A  large 
number  of  transport  aircraft  will  be  equipped  with  multiparameter  recording 
devices  to  collect  statistical  data  on  aircraft  loads  in  actual  service 
conditions . 

The  Netherlands  Civil  Aviation  Department  (RLD)  and  the  FAA  have  signed  a 
memorandum  of  cooperation  in  the  area  of  aircraft  structural  integrity,  with 
specific  reference  to  aging  aircraft.  As  part  of  this  cooperation,  the  National 
Aerospace  Laboratory  (NLR)  was  contracted  to  participate  in  the  flight  loads 
program.  One  of  the  tasks  specified  is  the  acquisition  and  review  of  European 
sources  of  flight  load  data  and  the  reduction  of  these  data  into  gust  statistics 
using  a  unified  reduction  procedure.  This  data  reduction  procedure  should  be 
compatible  with  the  analysis  format  to  be  applied  to  the  data  acquired  in  the 
current  FAA  flight  load  recording  program.  This  report  describes  the  procedures 
proposed  for  the  reduction  of  center  of  gravity  (vertical  acceleration  data’-) 
gust  statistics. 

In  the  development  of  these  procedures,  it  will  be  generally  assumed  that  a 
continuously  recorded  acceleration  time  trace  is  available  for  analysis, 
accompanied  by  additional  data  on  speed,  altitude,  weight,  aircraft 
configuration,  etc.  Specific  aspects  treated  successively  include: 

a.  The  separation  of  accelerations  due  to  gusts  and  maneuvers. 

b.  The  definition  of  peaks  and  valleys  in  the  acceleration  trace. 

c.  The  reduction  of  An  peaks/valleys  to  gust  velocities. 

It  will  be  proposed  to  simultaneously  apply  two  reduction  techniques:  one  based 
on  a  discrete  gust  model  using  the  so-called  Pratt  formula  and  the  other  based  on 
a  continuous  PSD-gust  model. 

Much  of  existing  old  load  data  bases  for  the  most  part  are  available  only  in  a 
very  reduced  format.  Methods  to  reevaluate  this  data  will  be  discussed  in  section 
6,  with  specific  reference  to  the  acquired  European  data  sources,  namely  NLR 
Boeing  747  ACMS  data,  ONERA-high  load  data,  and  RAE  fatiguemeter  data. 

The  report  concludes  with  a  summary  of  proposed  data  reduction  and  analysis 
procedures . 

2-t _ SEPARATION  OF  ACCELERATIONS  DUE  TO  MANEUVERS  AND  GUSTS. 

The  incremental  vertical  acceleration  at  the  center  of  gravity  (c.g.)  may  be  due 
to  either  maneuvers  or  gusts.  When  deriving  gust  statistics  from  acceleration 
traces  it  seems  logical  to  remove  the  maneuver- induced  acceleration  as  accurately 
as  possible  from  the  acceleration  history  and  consider  the  remaining  trace  as 
purely  gust-induced  acceleration.  In  the  past,  for  example,  in  reducing  the 
NACA/NASA  VGH  recordings,  this  separation  of  maneuvers  and  gusts  was  done 
manually  by  viewing  the  (analog)  record  traces.  Recently  however,  separation  was 


:  Appendix  A  gives  a  definition  of  "c.g.  vertical  acceleration"  as  meant 
in  the  present  report. 
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obtained  by  passing  the  acceleration  signal  through  a  variety  of  filters 
(reference  1)  . 

Figure  1,  reproduced  from  reference  1,  provides  an  example  of  the  type  of  result 
obtained  with  this  procedure,  illustrating  the  possible  ways  in  which 
accelerations  due  to  maneuvers  and  gusts  may  occur,  either  separately  or 
superimposed.  The  upper  figure  shows  the  unfiltered  trace.  The  second  figure 
shows  the  acceleration  signal  after  filtering  by  a  low-pass  filter;  the 
acceleration  trace  left  is  thought  to  be  due  to  maneuvering.  The  lower  figure 
shows  the  high-pass  filtered  signal,  showing  accelerations  due  to  gusts. 
Considering  this  figure  in  some  more  detail,  it  turns  out  that  the  complete 
separation  of  maneuvers  and  gusts  may  be  less  desirable  than  originally 
perceived . 

On  the  top  of  upper  figure,  five  typical  events  are  indicated,  namely  it 

1.  Shows  a  typical  acceleration  trace  due  to  a  turning  maneuver. 

-  The  induced  incremental  acceleration  is  always  positive,  and  has  a 
relatively  long  duration  (from  several  seconds  up  to  a  few  minutes) 

The  Typical  magnitude  is  in  the  order  of  0.15  g  (30°  bank  turn) . 

2.  Shows  the  same  type  of  turning  maneuver  with  gust  loads  superimposed. 
Typically,  one  may  think  of  turbulence  encountered  while  holding  at 
relatively  low  altitude. 

3 .  Shows  a  relatively  short  duration  acceleration  peak,  probably  due  to  a 
pitching  maneuver.  Positive  pitching  maneuvers  occur  typically  at  takeoff, 
during  rotation,  and  before  landing  in  flare.  The  magnitude  is  usually  limited 
to  the  order  of  0.10g.  Downward  pitching  maneuvers  occur  at  the  transition 
from  climb  to  cruise  and  from  cruise  to  descent.  The  magnitudes  of  these 
pitch-down  maneuvers  are  equally  small  or  probably  smaller  than  those  of  the 
pitch-up  maneuvers. 

4.  Shows  a  pure-gust  encounter.  The  acceleration  trace  has  a  highly  irregular 
character  with  relatively  high-frequency  content. 

5.  Shows  a  trace  where  maneuver  loads  are  superimposed  on  the  gust  loads. 

Note  that  the  maneuver  loads  are  both  upward  and  downward  and  of  a  rather 
irregular  nature.  Obviously  these  maneuvers  are  a  result  of  the  rurhu fence 
encounter;  the  pilot  applies  his  controls  to  counteract  the  effect  of  the 
turbulence.  The  latter  type  of  maneuvering  cannot  be  considered 
independently;  it  is  due  to  turbulence  and  should  be  considered  as  part  of  the 
aircraft  response  to  gust  input.  Hence,  one  may  decide  that  it  is  actually 
undesirable  to  filter  out  this  low-frequency  part  of  gust-induced  loading.  It 
should  be  noted  that  this  gust-induced  maneuvering  is  largely  in  pitch; 
corrective  turns  in  turbulence  remain  very  limited  indeed. 

The  majority  of  maneuver- induced  loads  associated  with  turning  maneuvers  can  be 
estimated  by  using  bank  angle  information.  If  the  bank  angle  is  recorded,  the 
acceleration  trace  may  be  corrected  for  turning  maneuvers  by  subtracting  a  An...r, 

«  (  l/cosip  -  1)  ,  where  !p  is  the  bank  angle  (reference  2)  . 

This  procedure  will  take  "'are  of  the  type  1  and  type  2  events  from  figure  1. 

As  explained  before,  the  majority  of  pitching  maneuvers  of  type  3  will  occur  at 
specific  points  in  the  flight  profile;  specifically,  the  maneuvers  at  rotation 
and  flare  can  easily  be  traced  and  removed  from  the  acceleration  trace  using  a 
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FIGURE  1.  FILTER  SEPARATION  OF  NORMAL  ACCELERATION  TIME  HISTORY 


suitable  detection  procedure.  To  summarize,  the  following  procedure  for 
separating  man-  vers  and  gusts  is  proposed: 

a.  Remove  the  acceleration  due  to  turning  maneuvers  by  subtracting  a  correction 
terr..  An.  ,,  =(l/cos$>  -  1)  from  the  acceleration  value. 

fc.  Remove  pitch  maneuver- induced  accelerations  at  specific  points  in  flight 
transition  such  as  takeof f /init ial  climb,  climb/cruise,  and 
approach/touchdown . 

c.  Do  not  remove  the  pitch  maneuver- induced  loads  in  turbulence,  as  they  are  part 
of  the  aircraft  system  response  to  turbulence. 

It  will  be  clear  when  reevaluating  existing  data  that  the  above  procedures  cannot 
or  can  only  in  part  be  applied  because  of  the  required  information,  e.g.,  the 
bank  angle  is  not  available.  The  only  advice  to  be  given  then  is  to  correct  as 
accurately  as  possible.  It  should  be  realized  that  in  any  case,  medium- rang  . 
transport  aircraft  turbulence  is  a  much  more  important  loading  source  than 
maneuvers  and  that  the  correction  for  maneuvers  is  not  of  vital  importance  for 
obtaining  valuable  statistical  data  on  gust  loads. 

3.  PEAK/VALLEY  SELECTION. 

Recall  that  the  objective  is  to  obtain  gust  statistical  data  from  recorded 
acceleration  traces  by  reducing  acceleration  peak/valley  values  to  derived  gust 
velocities.  The  first  step  necessary  in  the  procedure  is  the  selection  of  the 
peaks  and  valleys  to  be  reduced.  Figure  2  illustrates  the  procedure  generally 
applied  for  peak/valley  recognition  in  a  continuous  load  trace  to  avoid  the 
inclusion  of  peaks /val leys  associated  with  irrelevant  small  load  variations;  a 
range  filter  (R)  is  maintained  so  that  recognized  successive  peaks  and  valleys 
differ  at  least  R  in  value. 

If  we  wish  to  reduce  recognized  acceleration  peaks  to  upward  gusts  and 
acceleration  valleys  to  downward  gusts,  the  described  peak/valley  selection 
procedure  is  inadequate.  As  shown  in  figure  2,  acceleration  valleys  may  still 
have  a  positive  sign,  and  an  acceleration  peak  may  be  negative! 

For  this  reason  the  analysis  of  the  NACA  VGH  data,  a  more  drastic  peak/valley 
criterion  was  used.  This  is  shown  in  figure  3  and  usually  indicated  as  the  Peak 
between  Means  method.  Between  two  successive  crossings  of  the  mean  (in  case  of 
vertical  acceleration  the  1 . Og  -level)  only  one  peak  or  one  valley  is  recognized. 
To  further  ignore  irrelevant  load  variations  around  1 . Og  a  threshold  zone  of  An  = 
±  0.02  was  used  in  the  VGH  data  reduction,  as  shown  in  figure  4  (reference  3) . 

The  amount  of  peaks/valleys  that  are  rejected  by  the  Peak- between-Means  criterion 
depends  on  the  so-called  "irregularity"  of  the  signal.  The  irregularity  factor 
(k)  of  a  signal  is  defined  as  the  total  number  of  peaks  divided  by  the  total 
number  of  crossings  of  the  mean  in  the  positive  direction.  Figure  5a  shows  a 
signal  with  a  relatively  high  irregularity,  while  Figure  5b  depicts  a  typical 
narrow  band  signal  having  a  very  low  irregularity.  For  signals  with  an 
irregularity  factor  equal  to  1,  the  Peak-between-Means  criterion  does  not  reject 
any  intermediate  peaks  or  valleys,  simply  because  they  are  not  present  in  the 
signal . 

The  aircraft  response  to  turbulence  has  a  narrow-band  band  character  around  the 
pitch  response  frequency.  Figure  6  presents  a  part  of  a  c.g.  acceleration  trace 
recorded  in  a  Boeing  747  flying  through  turbulence  (reference  4),  showing  the 
narrow-band  character  and  associated  low  irregularity  of  the  signal.  Note  that 
application  of  the  Peak-between-Means  criterion  would  lead  to  recognition  of  all 
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FIGURE  2  CLASSIFICATION  OF  PEAKS  AND  TROUGHS  USING  A  RANGE  FILTER  R 


HGURE  3  PEAK-BETWEEN-MEANS  CLASSIFICATION  CRITERION 
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FIGURE  4.  METHOD  OF  EVALUATING  ACCELERATIONS  FROM  VGH  RECORDS, 
BY  THE  PEAK-BETWEEN-MEANS  CRITERION  (FROM  REFERENCE  3) 
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FIGURE  6.  EXAMPLE  OF  RECORDED  ACCELERATION  TRACE  (B747  IN 
TURBULENCE) 


FIGURE  7.  DISTRIBUTION  OF  RESPONSE  FOR  ONE  AND  TWO  DEGREES  OF 
FREEDOM  (FROM  REFERENCE  8) 


the  main  turning  points  in  the  load  sequence.  Hence,  the  Peak-between-Means 
criterion  is  less  drastic  than  may  appear  at  first  glance  when  applied  to 
aircraft  acceleration  data.  The  Peak-between-Means  criterion  is  well  suited  to 
be  applied  in  the  analysis  of  gust  -  induced  accelerations.  Also  in  order  to  remain 
compatible  with  the  existing  NACA  VGH  gust  data  base,  it  is  proposed  to  use  this 
criterion  in  the  reanalysis  of  existing  data  sources  as  well  as  in  future  data 
acquisition  projects. 


4.  REDUCTION  OF  An.- PEAKS /VALLEYS  TO  GUST  VELOCITIES. 

Two  methods  for  the  reduction  of  acceleration  peaks  to  derived  gust  velocities 
will  be  defined  based  on  essentially  different  representations  of  atmospheric 
turbulence.  The  first  is  based  on  a  discrete-gust  concept.  The  aircraft  c.g. 
response  to  gust  is  described  by  one  parameter,  indicated  as  C.  Recorded 
acceleration  peaks  are  reduced  to  derived  gust  velocities  Ud, .  The  gust  statistics 
that  will  eventually  be  obtained  describe  the  number  of  axceedings  per  kilometer 
of  Ude  for  a  number  of  altitude  bands. 

The  second  method  is  based  on  a  continuous-gust  concept.  The  aircraft  c.g. 
response  to  gust  is  described  by  two  parameters,  namely  A  and  N(0)  .  The  recorded 
acceleration  peaks  are  reduced  to  derived  gust  velocities  Ug. 

The  gust  statistics  obtained  initially  will  have  the  form  of  UQ-exceedance  curves 
for  a  number  of  altitude  bands.  These  curves  may  be  further  evaluated  to  yield 
for  the  various  altitude  bands  the  PSD-gust  parameters  P,,  Pz,  b,  and  b2 . 

It  is  proposed  to  use  both  methods  simultaneously  in  the  analysis  of  acceleration 
peak  data.  Thus,  two  sets  of  gust  statistics  will  be  obtained:  one  based  on  a 
classical  discrete-gust  model  and  easily  comparable  with  old  data  and  the  other 
set  based  on  a  more  advanced  and  realistic  gust  concept . 

a  - 1 , _ PI S .CRETE - SUS.T  APPROACH. ... 

The  atmospheric  turbulence  is  thought  to  consist  of  separate  discrete  bumps  with 
a  specific  gust  profile  and  magnitude  U.  The  discrete -gust  bump  induces  an 
aircraft  c.g.  acceleration  with  maximum  value  Anz  .  This  Anz  is  related  to  U  by: 

An., .  C  . U  (4.1) 

Equation  4.1  can  be  used  in  the  opposite  direction  to  calculate  the  magnitude  U 
of  the  gust  that  caused  a  recorded  incremental  load  factor  Anz : 

U  -  An Z/C  (4.2) 


The  response  function  C  is  generally  a  function  of  the  shape  and  length  of  the 
gust  and  the  response  properties  of  the  aircraft.  In  NACA  Report  1206  (reference 
5)  the  response  factor  C  was  calculated  under  the  following  assumptions: 

a.  The  gust  has  a  (1  -  cos)  shape  with  a  length  of  25-wing  chords. 

b.  The  aircraft  is  assumed  to  be  infinitely  stiff  and  to  respond  only  in 
plunge  (and  not  in  pitch) . 


Aerodynamic  inertia  is  included. 


The  result  found  for  C  could  be  approximated  by: 
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(4.3) 
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where 


F<V 


is  the  so-called  gust-alleviation  factor. 


Vi 


•68  Mg 

5.3  •  pg 


reading : 


(4.4) 


The  above  formula  is  the  well-known  so-called  "Pratt  formula",  derived  in 
reference  5.  It  has  been  used  extensively,  not  only  for  reducing  the  NACA  VGH 
data,  but  also  in  airworthiness  requirements  for  defining  gust  design  load 
factors.  The  physical  reality  of  the  Pratt  formula  is  rather  weak.  The 
assumption  of  a  gust  length  that  is  a  function  of  the  size  of  the  aircraft  flying 
is  inappropriate,  and  the  assumption  of  plunge  freedom  only  is  very  restrictive. 
The  great  advantage,  however,  of  the  Pratt  formula  is  its  simplicity. 

For  this  reason,  and  in  order  to  maintain  compatibility  with  the  existing  gust 
data  bases  which  have  been  obtained  using  the  Pratt  formula,  it  is  proposed  to 
use  this  formula  in  the  reevaluation  of  other  data  sources  and  in  the  analysis  of 
new  data . 


4.2  PSD  APPROACH. 

4 . 2-i  1 _ Detepmigatior} ,  of_JL- 

Rather  than  a  deterministic  sequence  of  discrete  bumps,  gust  activity  can  be 
described  as  a  random  process  of  atmospheric  turbulence.  The  c.g.  acceleration 
response  of  an  aircraft  flying  through  a  stationary  random  gust  field  with  rms 
value  ow  will  also  be  stationary  random,  with  rms  value  °An  The  ratio  0^n/0w  is 
called  A.  This  ratio  A  is  a  response  factor  that  is  to  a  large  extent  equivalent 
with  C  described  in  the  paragraph  4.1.  This  equivalence  has  led  to  the 
formulation  of  a  quasi -discrete  PSD-gust  velocity  U^,  for  which  holds 

Anz  .  A  .  Ua 


This  formulation  has  found  its  way  into  current  airworthiness  requirements,  where 
design  values  for  UQ  are  defined  as  a  function  of  altitude  (FAR  part  25,  Appendix 
G)  .  Full  discussion  of  PSD  methods  is  outside  the  scope  of  this  report  and 
reference  must  be  made  to  the  appropriate  textbooks.  At  this  point,  simply  recall 
that  A  is  found  by  integretation  over  all  frequencies  of  the  aircraft  transfer 
function  squared  times  the  normalized  power  spectral  density  function  of 
turbulence : 


A 


/  0“>!2 


(4 .5) 


Thus,  A  may  be  considered  as  a  weighted  average  response,  taking  into  account  the 
aircraft  response  characterist ics , that  is  the  sensitivity  for  specific 
wavelengths  (through  the  transfer  function  H)  on  the  one  hand  and  the  atmospheric 
properties  (which  wavelengths  do  occur?)  through  the  PSD  function  on  the  other. 
The  A  approach  is  based  on  a  much  more  realistic  turbulence  model.  Moreover, 
note  that  in  calculating  A  the  real  wavelength  distributions  as  occurring  in 
turbulence  were  considered  rather  than  assuming  an  aircraft  si ze - dependent  gust 
as  done  for  calculating  C. 
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A  may  be  calculated  using  the  same  simplifying  assumptions  with  regard  to 
aircraft  response  as  in  the  Pratt  formula,  that  is  infinite  stiffness  and 
response  freedom  in  plunge  only.  Such  calculations  have  been  made  in  the 
past,  for  example  in  reference  6  and  reference  7.  However,  as  pointed  out  in 
references  8  and  9,  the  assumption  of  no  response  in  pitch  is  unrealistic 
here.  The  weathercock  effect  associated  with  pitch  response  freedom  greatly 
reduces  the  aircraft  response  on  gusts  with  long  wavelengths.  Figure  7 
compares  the  aircraft  c.g.  acceleration  response  for  an  aircraft  with  plunge 
freedom  only  and  plunge  and  pitch.  (The  value  of  the  discrete-response 
function  C  is  not  very  dependent  on  the  inclusion  of  pitch  freedom  because  the 
assumed  gust  wavelength  of  25  chords  is  very  short) .  Hence,  a  simple  but 
realistic  expression  for  A  should  at  least  consider  plunge  and  pitch  response 
freedom.  Houbolt  (reference  8)  calculated  A  values  for  a  wide  variety  of 
aircraft, aircraft  conditions, and  flight  conditions.  Somewhat  to  his  surprise 
he  found  that  the  A  for  a  wide  range  of  aircraft  could  be  approximated 
accurately  by  using  the  following  expression: 


A 


Po  vs  CLa 

2  mg  /S 


F ( PSD  ) 


(4.6) 


where  F  (  PSD  ) 


11  .8 


\|  110 
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The  discrete-response 
compared  by  using  the 
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factor  C  and  the  continuous-response  factor  A  may  be 
expression  for  F(PSD)  with  F(pg)  as  given  in  (4.4): 


Note  that  in  FJPSD)  the  aircraft  response  is  described  by  two  parameters, 
namely  p,,  and  c.  The  parameter  c  is  used  in  combination  with  the  turbulence 
scale  L  and  defines  the  size  of  the  aircraft  in  relation  to  a  typical  gust 
wave  length.  Recall  that  in  the  Pratt  formula  such  a  parameter  is  avoided 
artificially  by  assuming  that  the  typical  wavelength  is  a  constant  times  the 
wing  chord  c. 

Figure  8  depicts  F(pq)  and  F(PSD)  as  a  function  of  pq  for  some  values  of  c/L. 
The  expression  given  above  for  A,  as  derived  by  Houbolt,  is  as  easy  to 
calculate  as  C  according  to  Tratt  and  has  in  our  opinion  a  much  better 
physical  background,  both  with  regard  to  the  turbulence  model  and  with  regard 
to  the  aircraft  response  model.  It  is  proposed  to  reduce  all  An  peaks  to 
Upvalues  using  the  above  Houbolt  formula. 

4.2.2  The  N ( 0 )  Problem. 

4 .2.2.1  Review  of  the  problem. 

According  to  the  discrete  gust  concept,  the  turbulence  consist  of  a  discrete 
number  of  gusts  or  bumps  per  kilometer;  each  aircraft  flying  through  this 
turbulence  encounters  the  same  number  of  bumps;  each  bump  Ud(,  causes  one 
load  bump  An  =  C.Ud  .  The  exceedance  curve  per  km  for  a  specific 
altitude/ flight  condition  has  the  form: 

An  An 

N  (An)  .  B,fl  "<7  •  B2e  V 
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Note  that  the  total  number  of  An  bumps  per  km  is  equal  to  B,  +  B 


For  the  PSD-gust  concept,  the  situation  is  different.  The  number  of  An  response 
peaks  depends  on  the  aircraft  response  properties  through  the  parameter  N(0) . 
Consequently,  the  exceedance  curve  has  the  form: 


N (An)  .  N(0) 


An  An 

Pje  *>>*  •  P2e 


Note  that  the  total  number  of  An  bumps  per  km  is  equal  to  N(Q) [P^  +  P^] ,  thus 
this  depends  on  the  aircraft  response  parameter  N(0) .  The  questions  to  be 
answered  are: 

a.  Should  the  N(0) -effect  be  included  in  the  An-data  reduction  procedures? 

b.  If  the  answer  is  yes,  how  should  this  be  done? 


N(0)  -  [per  km] 

nc 


or  N  [0)  -  No(0) 


[per  km ] 


where  Nq(0) 
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To  answer  the  first  question,  it  is  necessary  to  study  the  possible  variations  in 
N(0)  :  if  N{0)  would  be  approximately  equal  for  all  aircraft  types  and,  e.g.,  be 
only  a  function  of  altitude,  one  might  ignore  the  N(0)  effect,  or  rather  assume  a 
constant  N(0)  value,  for  example  in  NACA  4332.  (reference  6),  where  for  the 
conversion  of  discrete-gust  data  into  a  PSD  model,  fixed  N ( 0 )  values  of  6  per  km 
(10  per  mile)  and  4.7  per  km  (7.5  per  mile)  for  storm  and  non-storm  turbulence 
respectively  were  assumed.  To  study  this  N(0)  variation,  an  approximate  formula 
given  by  Houbolt  (reference  8)  for  N(0)  will  be  used  (see  Fig.  9) : 

Table  1  presents  calculated  No(0)  values  for  a  wide  range  of  commercial  aircraft, 
for  three  different  weight  configurations  namely  M.TOW,  operations  weight  empty 
and  an  average  weight.  It  may  be  noted  that  the  Nq(0)  values  cover  a  wide  range, 
from  20  for  the  smallest  aircraft  at  low  weight  to  about  4  for  the  largest 
aircraft  at  maximum  weight.  From  this,  the  variation  in  N(0)  appears  to  be  too 
large  to  be  ignored.  In  the  following,  a  method  to  include  N(0)  variation  in  the 
An-data  reduction  will  be  outlined. 


4 .2 .2. 2  Accounting  for  N(0)  variation 


The  following  method  is  proposed  to  include  N(0)  variation  when  reducing  An-  data 
to  UQ  data  (PSD-type  data) .  For  each  altitude  h,  a  reference  value  N(0)re,  is 
defined : 


N  ( 0) 


No(0) 
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F1GURE  9.  MU  ACCORDING  TO  HODBOLT  ,PSD  MODEL.  2-DOB 


wh -re  N  ihJ ...  is  a  fixed  value  and  p  is  the  density  pertaining  to  the  altitude  h. 
B.u'-'d  on  the  data  presented  in  table  1,  the  value  =  8  km  1  is  proposed 

(av.-rage  [or  long-  and  medium-haul  transport),  A  recorded  acceleration  peak  An 

is  reduced  as  follows: 

a.  thing  t:  he  values  of  V,h  and  aircraft  mass,  m,  at  the  instant  of  occurrence  of 
An,  the  response  parameters  A  and  N(0)  are  calculated  using  for  the 

cal cu 1  at:  ion  of  N{0)  the  previously  presented  Houbolt-formula . 

b.  One  recorded  An  peak  for  an  aircraft  with  N(0)  is  set  equivalent  to 

N ( O'  .  , 

- An  -  peaks 

,M  ( 0 ) 

tor  an  aircraft  with  N(0)re, . 

c.  The  value  An  is  reduced  to  gust  velocity  using 

"a  •  An /A 

In  nummary  one  recorded  An  occurrence  is  reduced  to 

N  (  0  ■  , . , 

- ,  ocourances  of  U„. 

N ;  v  ° 


In  this  way,  the  reduced  acceleration  data,  UQ,  from  different  aircraft  types  and 
different  weight  conditions  can  be  added.  The  resulting  statistic,  pertaining  to 
a  specific  altitude  band,  will  have  the  shape: 


Generally,  the  number  of  positive  gusts  larger  than  UQ,  N(U0+)  will  be  different 
from  the  number  of  negative  gusts  larger  than  U0,  N(U0~) . 

The  average  number  of  exceedings  of  U^,  N{UQ)  is  the  logarithmic  average  (see 
previous  figure) : 

N(|UJ)  -  *  n(uo) 


N  (U0)  (log  scale) 


5 .  CONVENTIONS 


In  order  to  establish  a  unified  procedure  for  reduction  of  acceleration  data 
resulting  in  gust  statistics  from  different  sources  that  can  easily  be  compared, 
it  is  desirable  to  establish  firm  conventions  on  a  couple  of  subjects. 

_ UNIT  SYSTEM. 

All  data  must  be  presented  in  Standard  International (SI)  units,  with  one 
exception,  namely  the  definition  of  altitude  bands  (see  section  5.3) .  The  main 
units  pertaining  to  the  SI  system  are: 


time 

second 

(s) 

length 

meter 

(m) 

mass 

kilogram 

(kg) 

force 

Newton 

(N) 

temperature 

degrees  Celsius 

CC) 

5^2 _ VELOCITY  DIMENSIONS. 

The  derived  gust  velocities  Ude  and  Ug  are  given  as  equivalent  velocities. 

As  a  consequence,  the  PSD-gust  parameters  b,  and  b2  will  also  be  given  in  terms  of 
equivalent  velocities. 

iJ _ SLIITODK.  BMPS,, 

The  derived  gust  statistics  (number  of  gusts  exceeding  Uds,  respectively,  U0,  per 
m  or  per  1000  m)  are  presented  per  altitude  band.  It  is  customary  in  aviation  to 
express  flying  altitudes  in  feet  rather  than  meters;  the  altitude  bands  will  be 
defined  in  feet  (as  an  exception,  see  section  5.1) . 

The  following  altitude  bands  will  be  used: 


< 

1500 

ft 

1500 

4500 

ft 

4500 

9500 

ft 

9500 

14500 

ft 

14500 

-  19500 

ft 

19500 

-  24500 

ft 

24500 

-  29500 

ft 

29500 

-  34500 

ft 

34500 

-  39500 

ft 

> 

39500 

ft 

_ TURBULENCE 

SCALE  L ■ 

In  accordance  with  current  airworthiness  regulations  (FAR  25,  JAR  25),  the 
turbulence  scale  L  is  assumed  to  be  762  m  (2500  ft)  at  all  altitudes. 


Note :  This  convention  is  usual  for  U„„,  but  the  existing  data  sources  often 
present  UQ  and  b,  and  bj  in  terms  of  true  velocity 
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fL _ REEVALUATION  OF  EXISTING  EUROPEAN  LOAD  DATA  SOURCES. 


In  the  previous  paragraphs,  a  general  procedure  for  reducing  c.g.  acceleration 
data  to  gust  load  statistics  has  been  established  on  the  assumption  that  a 
continuously  recorded  c.g.  acceleration  trace  were  available.  To  reevaluate 
existing  data  sources,  this  general  procedure  will  be  applicable  only  in  part 
due  to  the  fact  that  the  data  has  already  undergone  considerable  reduction.  In 
the  following,  the  format  of  available  European  Data  sources  will  be  reviewed  and 
methods  to  reduce  this  data  in  accordance  with  our  general  procedure  will  be 
defined . 

6.1  NLR  BOEING  B747  ACMS  DATA,  (reference  4) 

Acceleration  peak/valley  sequences  on  a  per  flight  basis  are  available  for  20,205 
flights.  Successive  peaks  and  valleys  have  been  determined  using  a  range  filter 
size  R=  0.18g.  The  following  information  is  available: 

1.  Flight  profile  data;  For  each  flight  the  speed,  altitude,  and  weight  as  a 
function  of  time  have  been  recorded. 

2.  Peak  data:  For  each  peak  or  valley,  the  Anz~value,  the  time  at  occurrence 
of  the  peak,  and  the  flap  setting  are  recorded.  For  a  limited  number  of 

flights,  the  bank  angle  at  peak-occurrence  is  also  recorded. 

This  implies  the  following  with  regard  to  the  applicability  of  our  reduction 
procedure : 

•  For  only  a  limited  number  of  flights,  the  acceleration  peak/valley  sequence 
can  be  corrected  for  (turning)  maneuvers  by  making  the  bank  angle  correction. 

•  The  (corrected)  peak/valley  sequence  may  be  further  reduced  using  the  Peak- 
between-Means  criterion. 

•  Sufficient  data  (speed,  altitude,  weight,  2nd  flap  setting  at  time  of 
peak/valley)  is  available  to  calculate  the  response  parameters  C,  A,  and  N(0) 
associated  with  each  peak/valley  so  that  gust  velocities  Ud>!  and  Ug  can  be 
derived . 

•  From  the  flight  profile  data  the  distances  flown  in  the  various  altitude 
bands  are  available  to  reduce  the  gust  exceedance  data  for  each  altitude  band 
to  gust  exceedings  per  km 

In  summary,  it  appears  that  the  general  data  reduction  procedure  is  applicable  to 
a  very  large  extent  to  the  B747  ACMS  data.  It  should  be  kept  in  mind  that,  as  the 
B747  ACMS  system  only  recorded  acceleration  traces  when  |An,  |  =  0 , 18g  was 
exceeded  (reference  4),  the  data  base  is  limited  to  gust  velocities  higher  than 
those  corresponding  with  this  incremental  acceleration  level  of  0.18g. 

6 . 2  ON ERA  HIGH  LOAD  DATA ■ 

The  ONERA  data  base  contains  information  about  relatively  high  loads,  namely 
acceleration  peaks/valleys  larger  than  |An,|  =  0.50g,  that  have  occurred  in  a 
very  large  number  of  flights  (about  800,000) .  These  peaks/valleys  have  been 
derived  from  ACMS- recorded  acceleration  traces  using  the  Peak-between-Means 
selection  criterion.  The  data  base  includes,  for  each  recognized  peak  or  valley, 
the  following  information: 
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1.  Incremental  load  factor  An*. 

2.  Aircraft  type,  tail  number,  and  date. 

3.  At  time  of  peak:  mass,  speed,  altitude,  flap  position,  and  C  . 

The  data  base  also  includes  the  total  number  of  flights  and  flight  hours  analyzed 
for  each  aircraft.  With  regard  to  applicability  of  the  standard  reduction 
procedures  we  may  conclude  that  sufficient  data  is  available  to  calculate  C,  A, 


and  N(G)  for  each  peak/valley  occurrance  and  hence,  to  calculate  exceedance 
curves  for  the  derived  gust  velocities  Ude  and  U^.  As  no  specific  flight  profile 
data  were  r-.-cot  ded,  estimated  average  mission  profiles  must  be  used  to  determine 
the  total  distances  flown  in  each  altitude  band  in  order  to  reduce  the  gust 
exceedance  data  to  exceedings  per  km.  Overall,  the  ONERA  data  appear  compatible 
with  our  standard  reduction  procedures. 

This  data  will  only  provide  information  about  severe  and  hence  relatively  rare 

au  s  t  s . 


<L .2  RAK  FATIGUE  METER  DATA.  (Reference  11) 

These  data  were  obtained  from  counting  accelerometers  measuring  center  of  gravity 
vertical  ,■»*••  1  •-■rat  ion .  The  instrument  counts  the  number  of  times  that  preset 

accelerat  '•  ci>  '.‘■  -‘■-Is  are  exceeded.  To  avoid  spurious  counts  due  to  small 
acc-eler.V  :  •  >n  variations,  the  exceedance  of  a  level  is  only  counted  after  the 
acceleration  has  dropped  below  another  predetermined,  lower  level.  More  details 
on  the  counting  procedure,  in  comparison  with  the  Peak-between-Means  method,  are 
described  it.  appendix  C. 

At  .  istorv  mi  of  a  few  minutes  the  counters  were  photographed  together  with  an 
altimeter,  speed  indicator,  and  clock.  The  data  available  cover  about  26,000 
hour:.:  of  flvirig  and  include  a  variety  of  aircraft.  The  current  data  base  is 
available  on  magnetic  tape  and  is  obviously  the  result  of  considerable  data 
proc-visi ng,  starting  from  the  original  information  contained  on  photographic 
images.  Available  data  are  presented  in  the  following  format: 

a .  For  each  aircraft  type,  a  large  number  of  intervals  are  distinguished  with 

regard  to: 

•  Flight  C'-Vi  -nt  (e.g.,  climb,  cruise,  or  descent) 

•  Altitude  (a  set  of  altitude  bands  is  defined) 

•  Speed  (a  number  of  speed  ranges  are  defined) 

•  We  I  glit  (a  number  of  weight  ranges  are  defined) 

b.  For  each  interval  (one  flight  segment,  altitude  band,  speed,  and  weight 
range),  hhe  following  information  is  presented: 

•  The  total  time  spent  and  distance  flown  within  the  interval. 

•  The  total  number  of  crossings  of  the  various  acceleration  counting  levels. 

In  a .tier  to  make  these  data  compatible  with  the  standard  reduction  procedure,  the 
count' vd  b  v-o  crossings  should  be  converted  to  counted  peaks  and  valleys.  As 
town  in  appendix  C,  the  difference  in  number  of  crossings  of  two  successive 
(po-'ot  i  ■/•:•)  levels  is  equal  to  the  number  of  peaks  with  magnitude  larger  than  the 
Iow-t  and  orroiiler  of  the  higher  of  these  two  levels.  In  appendix  B  discrete 
e-pi  i  v  «•)  *-nt  peak  values  for  all  peaks  between  two  level  cross  counting  levels  have 
he<  n  q.,f  j  nr-d .  Using  these  values,  the  counted  level  crossings  for  each  interval 
'•■•in  rKw  bo  converted  to  counted  peaks/valleys  for  each  interval. 
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Information  (speed,  altitude,  and  weight)  is  available  to  calculate  C,  A,  and 
N(0)  and  hence,  to  reduce  the  acceleration  peaks/valleys  to  derived  gust 
velocities  U,..  and  U^.  The  speed,  weight,  and  altitude  are  only  known  within  a 
relatively  wide  range;  hence,  the  relative  accuracy  of  the  derived  gust 
velocities  will  be  lower  than  for  the  previously  described  NLR  and  ONERA  data . 

In  addition,  accelerations  due  to  maneuvers  cannot  be  eliminated  and  that,  as 
explained  in  appendix  C,  the  applied  counting  technique  may  lead  to  a  slight 
overestimation  of  relevant  peaks. 

In  summary,  we  may  conclude  that  the  Fatiguemeter  data  can  be  reduced  and 
presented  such  that  it  is  compatible  with  the  other  data  sources. 

7,  SUMMARY  OF  PROPOSED  REDUCTION  PROCEDURES. 

A  summary  is  presented  of  the  successive  steps  in  the  general  procedure  to  reduce 
a  recorded  acceleration  trace  to  gust  velocity  data. 

7  ■  l _ Sorting  points 

We  start  with  a  continuous  record  of  vertical  acceleration  at  or  near  the 
aircraft  center  of  gravity.  It  is  assumed  that  the  signal  is  adequately 
filtered  to  eliminate  structural  responses.  If  the  signal  is  available  in  digital 
form,  the  signal  should  in  any  case  be  free  of  frequency  components  higher  than 
half  the  sampling  rate. 

.7. -.2 . Ell.mipa.Li.9n,,.  o£.,,m.a.0.e.UVeX.-I.OAds.  ■ 

If  the  bank  angle  history  is  available,  the  acceleration  history  should  be 
corrected  for  turning  maneuver- induced  load  factor  increments  by  subtraction  of  a 
correction  term:  An,.,,  =  (l/cosq)-l)  .  If  specific  pitch  maneuvers  can  be 
recognized  in  the  signal  (e.g.,  rotation  at  takeoff,  end  of  climb,  etc.)  the 
induced  load  factor  increments  should  be  removed  from  the  acceleration  trace. 

Z_i3 _ Peak/vallev  selection. 

The  corrected  (i.e.,  maneuver  removed)  acceleration  history  is  searched  for  peaks 
and  valleys  using  the  Peak- be tween  -  Means  selection  criterion. 

7.4  Reduction  of  accelerations  to  gust  velocities. 

The  classified  acceleration  peaks  and  valleys  An,  are  reduced  to  derived  gust 
velocities,  following  a  discrete-gust  approach  and  a  PSD-gust  approach 

respectively . 

7.4.1  Discrete  Gust:. 

Each  single  An,  is  reduced  to  one  gust  U according  to: 


C 


with  :  C 


4  mg  /  £ 


(7.2) 


.1  9 


(7.3! 


where 


•88m, 


5 . 3 


7.4.2.  PSD  Gust 


single  An,  is  reduced  to 


I  j 

o 


No(0i ref 

- - - /gusts  /with /magnitude  /UQ 


according  to: 


(7.4. 


with  A 


PqVeC^ 

2  mg  /S 


-F ( PSD  ! 


(7.5) 


where  F ( PSD  ) 


11.8  f  5  U 

/n  {  2lJ  110  .  Mg 


(7.6) 


and 


N  (Oi 


lie 
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SPo  Ci*  c 


(7.7) 


l,t>.  •'■•it  EXCEEDANCE  CURVES  . 

The  following  procedure  applies  to  both  Ude  and  U0  data,  but  will  be  described 
here  for  U,„  data  only. 

7 . f . l  Altitude  Bands. 

The  gust  peaks/valleys  are  sorted  according  to  the  altitude  band  in  which  they 
occur.  The  altitude  bands  have  been  defined  in  chapter  5.3. 

In  addition,  the  total  distance  flown  within  each  altitude  band  is  determined. 

7,5.2  Gust  Exceedance  Data, 

For  each  altitude  band  peak/valley, data  areacumulated,  resulting  in: 

•  For  a  set  of  positive  values:  -  number  of  gust  peaks  larger  than  that 

value 

•  For  negative  U,,„  values:  -  number  of  "gust  valleys"  lower  than  that  value 

Results  are  presented  as  total  numbers  and,  after  dividing  by  the  distance  flown 
within  the  altitude  band,  as  number  per  kilometer. 
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7-6  PSD-GUST  PARAMETERS 


Determine  the  average  number  of  exceedings  (up  or  down)  of  U0  from 
N(|Ue|)  -  y/wfj 4  • 


Determine  for  each  altitude  band  the  parameters  Plt  P2,  b, ,  and  b2  by  fitting  the 
curve 


-  N<°)  , 


P2e 


■e; 


to  the  derived  N(|U_|)  curve. 


Here  ,  N  (0)  ref  -  8  x 


■46 

-^1  km' 


where  p  is  the  air  density  pertaining  to  the  midvalue  hml(1  of  the  altitude  band 
under  consideration. 


7.7  REEVALUATION  OF  EXISTING  GUST  LOAD  DATA.  SOURCES ■ 

When  reevaluating  existing  gust  load  data  sources,  it  is  important  to  make  these 
data,  as  much  as  possible,  compatible  with  those  obtained  using  the  above 
reduction  procedures.  Key  elements  are  Peak-between-Means (peak/valley 
recognition) and 

expressions  used  for  C,  A  and  N(0)  . 

The  counting  accelerometer  data  are  yield  level  cross  data.  These  should  first 
be  reduced  to  equivalent  peak/valley  data  using  the  principles  outlined  in 
appendix  B. 

8 .  CONCLUSIONS . 

•  A  general  procedure  has  been  defined  for  the  reduction  of  measured  c.g. 
accelerations  to  derived  gust  velocities, 

•  Accelerations  due  to  turning  maneuvers  may  be  eliminated  if  the  bank  angle  ...» 
been  recorded.  In  addition,  the  accelerations  due  to  specific  pitch  maneuvers 
can  be  eliminated  if  they  can  be  recognized  in  the  flight  trace. 

•  If  the  acceleration  data  are  reduced  on  the  basis  of  a  continuous-gust 
concept,  the  variations  of  the  response  parameter  N(0)  should  be  included. 

A  method  to  do  so  has  been  defined. 

•  Available  European  gust  data  sources  can  be  reevaluated  and  reduced  into  a 
format  compatible  with  the  defined  analysis  procedures. 
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APPENDIX  A 

DEFINITION  OF  C.G.  VERTICAL  ACCELERATION 

The  equations  of  motion  of  an  aircraft  may  be  expressed  with  reference  to  an 
aircraft  fixed  stability  axis  system  (x,y,z) .  The  origin  of  this  axis  system  is 
the  aircraft  center  of  gravity  and  the  x  axis  is  defined  in  such  a  way  that  at 
time  t=o  the  x  axis  coincides  with  the  direction  of  the  aircraft  velocity. 

In  case  of  symmetric  motion,  the  equations  of  motion  of  an  elastic  aircraft  can 
be  described  as  follows: 

m(w-U.q)  .  Z  (A.  i; 


Iy  q  -  L 


(A.  2) 


where  U 


w 

q 

m 

iy 

and  ; 


velocity  of  c.g.  in  x  direction 
velocity  of  c.g.  in  z  direction 
rotational  velocity  around  y  axis 
aircraft  mass 

moment  of  inertia  around  y  axis 


HuU(t)  ♦  KuU(t>  -  Q4  (i-3, 


(A. 3) 


with  M,, 

Ki. 

Ut) 

Qi 


=  generalized  mass  i'h  mode 
=  generalized  stiffness  ith  mode 
=  generalized  displacement  ith  mode 
=  generalized  force  ict>  mode 


Equations  (A.l)  and  (A. 2)  describe  the  motion  of  the  aircraft  center  of  gravity. 
The  (m~2)  equations  (A. 3)  describe  the  deformation  of  the  structure  under  loading 
with  respect  to  the  aircraft  axis  system,  (x,y,z).  The  displacement  function  of 
the  i'-1'-  elastic  mode  reads 

^i(X/y>z»t)  -  (x,y,zKj(t) 


where  <&,  is  a  vector  with  components  and  in  x,y  and  z  direction 

respectively.  The  displacement  of  a  point  (x'.y'.z1)  in  the  structure  in  z 
direction  at  time  t,  W,  (x‘ , y' ,  z1 , t )  is  equal  to 

HjlxSy'.z'.t)  *  £  <frzl(x  \  y  \  z  Mf^t) 


The  c.g.  vertical  acceleration  az  is  defined  equation  A.l 

a  z  -  w  -  U  .  q 


(A.  5) 


This  is  the  quantity  we  would  like  to  measure  with  an  on  board  accelerometer. 
Unfortunately,  the  onboard  accelerometer  is  fixed  to  the  (deforming)  structure 
and  does  not  remain  in  the  c.g.,  which  is  in  fact  a  virtual  point;  the 
acceleration  actually  measured  by  the  accelerometer  mounted  in  point  x=0,  y=0. 


A- 1 


2=0,  is  equal  to 


a  j  -  w  -  U-q  .  <t2l  ( x.y.z-0 )  if 1 ( 1 ) 


iu 


(A.  6) 


Figure  A-l  presents  calculated  weighted  power  spectra  (w)  for  a  short-haul 
jet  transport  flying  through  stationary  turbulence.  The  spectra  refer  to  the 
real  c.g.  vertical  load  factor  increment  An,  =  a,/g  and  the  measured  value  An*  = 
a  */  9 

The  spectra  for  An.,  have  been  calculated  for  three  cases,  namely: 

1.  The  aircraft  responds  only  in  plunge  (1  DOF) 

2.  The  aircraft  responds  in  pitch  and  plunge,  but  is  assumed  to  be  infinitely 
stiff  (2  DOF) 

3.  As  in  2,  but  the  aircraft  flexibility  is  taken  into  account  by  8  flexible 
modes  (10  DOF) 


We  note  that  inclusion  of  flexibility  has  a  certain  but  limited  effect  on  the  An, 

response:  the  response  peak  shifts  slightly  to  the  lef t ( frequency) ,  and  the 
overall  response’  is  slightly  lower.  Note  also  that  the  response  is  largely 

restricted  to  say  18  rad/s  or  3  Hz;  the  energy  contained  in  higher  frequencies  is 
negligibly  small.  Looking  now  at  the  spectrum  pertaining  to  the  measured  load 
factor  increment  Anz,  we  note  that  the  spectrum  starts  to  deviate  from  that  of 
An,  above  that  frequency.  The  Anz~ spectrum  contains  peaks  at  high  frequencies, 
corresponding  with  structural  response  frequencies. 


♦ 

In  order  to  obtain  the  correct  An,  signal,  the  "measured"  Anz  signal  must  be 
filtered  so  that  the  high-frequency  components,  above  3  Hz,  are  eliminated. 


The  sample  frequency  of  8  Hz  often  used  in  recording  c.g.  vertical  accelerations 
in  transport  aircraft  is  reasonably  compatible  with  this  cut-off  frequency  of  3 
Hz.  (Sample  frequency  must  be  at  least  two  times  the  highest  frequency  present, 
to  avoid  fold-over  effects.) 


* 


2 

Note:  We  may  recall  that  the  overall  Anz  response  az,  is  proportional  to  the 
area  under  the  curve  drawn: 

m  m 

oL-  f*^  (w)dm.  J  d(  log  a) 

o  **o 


A-  2 


t - T - f — i — l"  'l"  rr-y  .  i - ! - 1 — t — r 


' - plunge  only 

- PrTCH  ♦PLUNGE 

_■ -  PITCH,  PLUNGE  ♦»  EM 

- PITCH.  PLUNGE  EM 

WEIGHTED 

POWERSPECTRUM 

M 


AIRCRAFT :  m  *  40000  teg 
H  •  7000  m 
V  *  230  m/i 


100 


w  (rad/t) 

FIGURE  Al.  WEIGHTED  POWERSPECTRUM  OF  DUE  TO  TURBULENCE  FOR  TWIN- 
ENGINED  SHORT  HAUL  TRANSPORT  AIRCRAFT 
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APPENDIX  B 

REPLACEMENT  OF  A  LEVEL  CROSSING  BY  AN  EQUIVALENT  PEAK  VALUE 


A  counting  accelerometer  presents  the  number  of  times  that  specific 
acceleration  levels  An,  are  exceeded.  Suppose  that  the  number  of  exceedmgs 
of  level  An,  is  one  higher  than  the  number  of  exceedings  of  the  next  higher 
counting  level  An,,,  =  An,  +  dAn.  This  implies  that  one  peak  with  magnitude 
between  An_,s  and  An,  +  dAn  has  occurred.  In  the  following  a  discrete- 
equivalent  value  for  the  positions  of  this  peak  will  be  derived. 


The  exceedance  curve 
for  An,  over  a  limited 
An  range,  can  be  fairly 
well  approximated  by  a 
straight  line  in  a  semi- 
logarithmic  grid. 

N(An)  =  ce  3 (B.  1) 


The  total  number  of  peaks  between  An,  and  An,„  =  An,  +  dAn,  is  equal  to 
N(An,)-  NfAn^).  [l  -  e . 


( B  .  2  ) 


The  probability  that  a  peak  within  this  interval  is  smaller  than  An,  +  dAn*  is 
given  by: 


P  ( dAn*) 


^  ^  S .  dAn* 
1  -  e  -s-'S" 


0<An*<An 


(B.  3} 


We  will  derive  equivalent  peak  values  An,  +  dAn„  on  the  basis  of  two  different 
criteria : 

1.  The  mean  or  most  probable  peak  value;  dAn„  =  EfdAn*] 

2.  dAne  is  the  "median"  value  dAn,  that  means  P(dAnJ  =  H. 

Ad  1: 

An  * 

Rewriting  (B.3)  with  B  =  S . An  and  x  =  yields. 

1  -p,  Bx 

P(x)  -  —  0<x<  1  (B.4) 

1  -  e  'R 


p(x) 


dP  Be  814 
dx  1  -  e  “B 


Osxsl 


x<0  and  x>l 


B-l 


«  1 

U(x)  E  [  x  ]  -  j  p(x)-x-dx  -  j 

-*>  O 


Be 

(1  -  e'B) 


dx  - 


1  -  e  'B(l  .  B) 
B(l  -  e  B) 


(B.  5} 


Ad  2: 


1  - 


-  e 


e 
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.  e 


x  - 


B 


in 


(B.  6) 


Numerical  values: 


The  mean  p(x)  and  median  x  according  to  (B.5)  and  (B.6)  respectively  are  a 
function  of  B:  B  =  S.dAn  depends  on  the  slope  S  of  the  An  exceedance  curve 
and  the  difference  dAn  between  successive  counting  levels.  A  review  of 
recorded  acceleration  data,  presented  in  references  4,  10  and  11  revealed 
that  the  value  of  S  ranges  from  about  6  for  a  steep  An  spectrum  to  about  15 
for  a  flat  spectrum.  For  example,  figure  B-l  shows  recorded  acceleration 
spectra  for  the  F-28  aircraft,  yielding  an  average  S  value  of  about  8.8. 

The  distance  between  successive  counting  levels  used  in  the  RAE  fatiguemeters 
(Reference  section  1.2)  is  either  dAn  =  0.1  or  dAn  =  0.2.  Thus,  the  total 
range  for  B  values  is  from  0.1  *  6  =  0.6  to  0.2  x  15_  =  3.0.  Table  B-l  gives 
values  for  the  mean  p(x)  (e.g.,  B.5)  and  the  median  x  (e.g.,  B.6)  as  a 
function  of  B.  For  all  B  values,  the  difference  between  mean  and  median  is 
relatively  small.  Over  the  B  range  considered,  the  mean  decreases  from  0.458 
to  0.28  and  the  median  from  0.438  to  0.21.  On  the  basis  of  these  findings, 
it  has  been  decided  to  assume  average  equivalent  x  values  of  0.40  and  0.33 
for  An  =  0.10  and  An  =  0.20  respectively. 

Table  B-2  presents  the  resulting  equivalent  peak  values  associated  with  the 
crossing  levels  for  the  two  types  of  fatigue  meters  used  in  the  RAE 
measurements  of  reference  11. 
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FIGURE  B1  EXAMPLE  OF  RECORDED  ACCELERATION  SPECTRA 
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TABLE  B-l 


MEAN  (i(x)  AND  MEDIAN  x  AS  A  FUNCTION  OF  B  =  S*dAn 


B 

mean 

aadian 

0 . 5 

0.45851 

0.43814 

0 . 6 

0.45030 

0.42610 

0 . 7 

0.44214 

0.41423 

0 . 8 

0 .43403 

0 .40256 

0 . 9 

0.42599 

0.39110 

1 . 0 

0.41802 

0 .37989 

1 . 1 

0.41013 

0.36892 

1.2 

0.40232 

0.35822 

1.3 

0.39460 

0.34780 

1.4 

0.38697 

0.33766 

1 . 5 

0 . 37945 

0.32782 

1 . 6 

0.37203 

0.31828 

1.7 

0.36472 

0.30904  | 

1 . 8 

0.35752 

0.30009 

1 . 9 

0.35044 

0.29145 

2 . 0 

0.34348 

0.28311 

2 . 1 

0.33665 

0.27506 

2.2 

0 . 32994 

0 . 26730 

2  .3 

0.32335 

0 .25983 

2 .4 

0.31690 

0.25263 

2 . 5 

0.31057 

0.24570 

2 . 6 

0.30438 

0.23904 

2 . 7 

0.29832 

0.23263 

2 . 8 

0 .29240 

0.22647 

2 . 9 

0 . 28660 

0.22055 

3 . 0 

0 .28094 

0.21485 

TABLE  B -  2 


EQUIVALENT  PEAK  LEVELS  TO  BE  USED  IN  THE  REEVALUATION 
OF  RAE  FATIGUEMETER  DATA 


MECHANICAL 

INSTRUMENTS 

ELECTRICAL 

INSTRUMENTS 

cross  level* 

equivalent  peak 

cross  level* 

equivalent  peak 

An 

(An+cLAn) 

An 

(An+cLAn) 

0 .23 

0.27 

0 .20 

0.24 

0 .33 

0.37 

0.30 

0 . 34 

0.43 

0.47 

0.40 

0.466 

0 . 52 

0 . 56 

0.60 

0.666 

0 . 62 

0.66 

0.80 

0 .866 

0.72 

0 . 76 

1.00 

1 . 066 

0 . 82 

0.86 

1.20 

1 . 266 

0 . 92 

0 . 96 

1.40 

1.466 

1 . 02 

1.10  1,1 

1 .60 

1.68  m 

Notes:  *:  from  reference  11,  table  1 

(1) :  more  or  less  arbitrarily  chosen 
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APPENDIX  C 

THE  COUNTING  PRINCIPLE  OF  THE  FATIGUEMETER 

Peak  count  methods  and  level  cross  count  methods  are  closely  related.  The 
difference  in  number  of  crossings  of  two  successive  levels  is  equal  to  the 
number  of  peaks  minus  the  number  of  valleys  between  these  two  levels 
(reference  12).  If  the  signal  in  question  has  a  narrow-band  character,  as  in 
the  case  of  c.g.  accelerations  due  to  turbulence,  the  number  of  valleys 
between  two  positive  levels  and  the  number  of  peaks  between  two  negative 
levels  will  approach  zero,  and  the  difference  in  crossings  of  the  two  levels 
will  be  equal  to  the  number  of  peaks  if  the  two  levels  are  positive  and  equal 
to  the  number  of  valleys  if  the  two  levels  are  negative.  Figure  C-l 
*  illustrates  the  counting  principle  of  the  fatigue  meter.  A  counting  level  is 

cocked  when  the  signal  crosses  that  level.  The  counting  of  that  level 
crossing  is  completed  when  the  signal  passes  another  reset  level,  which  is 
■»  closer  to  lg.  This  counting  principle  is  documented  in  the  literature  as 

restricted  level  cross  counting.  If  the  reset  levels  for  all  counting  levels 
were  set  at  lg  (An  =  0) ,  the  Peak-between-Mean  result  could  be  directly 
obtained  from  the  counted  level  crossings:  the  number  of  peaks  between  level 
j  and  level  j+1  would  be  equal  to  the  difference  in  crossings  of  the  two 
levels,  even  if  the  signal  does  not  exhibit  a  narrow-band  character. 

Table  C-l,  from  reference  11,  shows  that  the  reset  levels  maintained  do  not 
fully  comply  with  this  criterion.  This  means  that  more  peaks  and  valleys 
might  be  recognized  than  according  to  the  Peak-between-Means  criterion.  From 
the  example  of  figure  C-l,  two  crossings  of  level  j+1  and  one  of  level  j+2 
are  counted,  to  be  interpreted  as  one  peak  above  level  j+2  and  one  between 
j+1  and  j+2.  The  Peak-  between-Means  criterion  would  only  recognize  the  one 
peak  above  level  j+2  and  ignore  the  one  between  j+1  and  j+2.  Since  the 
acceleration  signal  is  of  narrow  band  nature  and  the  reset  levels  for  the 
most  relevant  counting  levels  are  very  near  the  mean (An  =0),  it  is  felt  that 
this  effect  may  be  ignored  and  that  differences  in  counted  crossings  of 
successive  positive  levels  may  be  interpreted  as  number  of  peaks  between 
these  levels  (and  valleys  in  case  of  successive  negative  levels.) 
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TABLE  C-l 

ACCELERATION  INCREMENT  LEVELS  (G)  COUNTED  BY  THE 
RAE  FATIGUE  METERS  (ELECTRICAx.  INSTRUMENTS) 


COUNTER  COCKED 

COUNT  COMPLETED 

0.2 

0.0 

0.3 

0.0 

0.4 

0.1 

0.6 

0.2 

0.8 

0.3 

1 . 0 

0.4 

1.2 

0.6 

1.4 

0.8 

1.6 

1 . 0 

from  reference  11 
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